Two-dimensional frequency-domain visco-elastic full waveform inversion: Parallel algorithms, optimization and performance

Full waveform inversion (FWI) is an appealing seismic data-fitting procedure for the derivation of highresolution quantitative models of the subsurface at various scales. Full modelling and inversion of viscoelastic waves from multiple seismic sources allow for the recovering of different physical parameters, although they remain computationally challenging tasks. An efficient massively parallel, frequencydomain FWI algorithm is implemented here on large-scale distributed-memory platforms for imaging two-dimensional visco-elastic media. The resolution of the elastodynamic equations, as the forward problem of the inversion, is performed in the frequency domain on unstructured triangularmeshes, using a low-order finite element discontinuous Galerkin method. The linear system resulting from discretization of the forwardproblem is solvedwith a parallel direct solver. The inverse problem,which is presented as a non-linear local optimization problem, is solved in parallel with a quasi-Newton method, and this allows for reliable estimation ofmultiple classes of visco-elastic parameters. Two levels of parallelism are implemented in the algorithm, based onmessage passing interfaces andmulti-threading, for optimal use of computational time and the core-memory resources available on modern distributed-memory multicore computational platforms. The algorithm allows for imaging of realistic targets at various scales, ranging from near-surface geotechnic applications to crustal-scale exploration. & 2010 Elsevier Ltd. All rights reserved.